Measuring non-linear functionals of quantum harmonic oscillator states. 
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Using only linear interactions and a local parity measurement we show how entanglement can be 
detected between two harmonic oscillators. The scheme generalizes to measure both linear and non- 
linear functionals of an arbitrary oscillator state. This leads to many applications including purity 
tests, eigenvalue estimation, entropy and distance measures - all without the need for non-linear 
interactions or complete state reconstruction. Remarkably, experimental realization of the proposed 
scheme is already within the reach of current technology with linear optics. 



In the context of quantum communication and com- 
puting protocols, measures such as purity, fidelity and 
entanglement characterize the performance and non- 
classical resources in a physical experiment. Such mea- 
sures not only provide a link with theoretical models, 
but also provide standards for defining benchmarks 0. 
With continuing technical developments in these areas, 
it is becoming increasingly necessary to identify practical 
and efficient schemes to measure such quantities. 

One obvious method is to first reconstruct the com- 
plete density matrix from a series of measurements us- 
ing, for example, the well-know procedure of quantum 
state tomography From this reconstructed den- 

sity matrix the desired measure can then be computed. 
Although this technique is realizable, it is not efficient 
as much more information about the quantum state is 
obtained than is actually needed. 

A more direct method was recently proposed by Filip 
H , and expanded upon by others [a, 0, Q • They pro- 
posed a specific quantum circuit capable of measuring 
the non-linear functional Tr(p^p^) of the density matri- 
ces p a and p b . The scheme requires as inputs, k and 
I copies of the states p a and ph, respectively. It was 
independently shown by Brun [2j that any polynomial 
function of a state up to degree q can be estimated by a 
joint measurement on q copies of the system. Amazingly, 
when the systems are entangled the measurement dou- 
bles as an entanglement witness 0, giving a negative 
value for some entangled states while ensuring a positive 
value for all separable states. 

The circuit, although elegant, unfortunately requires a 
nontrivial interaction between a control qubit and the k 
and I copies of each system - the targets. In its simplest 
form with two targets, k + I = 2, the interaction must 
generate a controlled-SWAP operation between the con- 
trol and the two targets. Such an operation does nothing 
to the targets if the control is in the logical zero state 
and applies the unitary swap operation 



v\<t>i)\<h) = \<h)\4n 



(i) 



to the two targets if the control is in the logical one 
state, where here |</>i) and \(j> 2 ) are arbitrary target states. 
Implementing this gate is a challenging requirement for 



any experimental architecture, especially when the tar- 
get states are of arbitrary dimension. In the case of a 
harmonic oscillator where the dimension of the Hilbcrt 
space is infinite, implementing this gate usin g, fo r exam- 
ple, non-deterministic KLM-type gates [ill Il2l] is only 
possible when there is at most one quanta of energy in 
each oscillator. 

In this paper we propose an alternative to the Filip 
scheme for harmonic oscillators which, remarkably, re- 
quires only linear coupling between different oscillators 
making it realizable with current technology. Like the 
Filip proposal, the scheme can measure both linear and 
non-linear functionals of the density matrix, while for 
entangled oscillators the scheme acts as an entanglement 
witness. We begin by first reviewing some of the prop- 
erties of the swap operator @, 0, H, @ , along with its 
generalization to multi-particles, and show how it relates 
to a host of relevant information measures. 

For the swap operator defined in Eqn. 0) it is straight- 
forward to show that for two separable states p a and p , 

Tr(p aPb ) = Tr ab (p a ®p b V). (2) 

Although this seems like a trivial relation, it suggests 
that a direct way to measure the overlap between two 
unknown states is to measure the expectation value of the 
swap operator (V). In the case where both systems are 
in the same state pi, the expectation value is equivalent 
to the purity of the system, (V) = Tr(pf). From the 
purity and the overlap one could obtain, for example, 
the Hilbert-Schmidt distance Tr[(p a — Pb) 2 ] ■ 

In the case where the state is not separable Eqn. © is 
no longer valid. To illustrate this take, for example, the 
entangled state (|a)|/3) — \{3)\a))/y/2, where \a) and \ f$) 
are orthonormal. From 0) the expectation value of the 
swap operator is —1 which, in contrast to Q, is negative. 
This is a example of the se par ability criteria of Horodecki, 
Horodecki and Horodecki 13] which states that a density 
matrix p is entangled iff there exists a Hermitian operator 
H, an entanglement witness, such that 

Tr(pH) < 0, (3) 

while for all separable states p sep 

Tr(PsepH) > 0. (4) 



From Eqn. J5J and by explicit example we see that the 
swap operator is an example of an entanglement witness. 
A measurement of (V) is said to have witnessed the en- 
tanglement when the outcome (V) < 0. In fact, it can 
be shown that the swap operator is an optimal entan- 
glement witness in the sense that it forms a hyper-plane 
that is tangential to the convex set of separable states 
[T3, Hij | . That is a separable state p' sep exist such that 
Ti[p' sep V] = 0. Such a state is \a) <g> |/3). 

A generalization of the swap operator to multiple sys- 
tems is defined as 

v N \(h)\4>») ■ ■ ■ \<Pn) = ■ ■ ■ \4>n-i), (5) 

which is not Hermitian for N > 3. For separable states 
Psep = Pi ® P2 <S> ■ • ■ Pn Eqn. (J2J generalizes to 

Tr[p sep V N ] = Tr[p a p b ... p N }. (6) 

For N identical copies of a state, abbreviated as p® N , 
this becomes 

Tr[ P ® N V N ]=Tr(p N ) = J2^ (7) 

i 

where Xi is the i th eigenvalue of p. We note that for a set 
of N identical pure states the expectation value is unity, 
while if the set is not homogeneous the expectation value 
will be less than unity. This serves as a practical state 
discrimination test of multiple systems |23j . 

For a rf— dimensional system the spectrum of p can 
be obtained from the (d — 1) values of Tr(p N ) for N = 
2,3, ... ,d JQj. Once the spectrum is known any nonlin- 
ear functional of the general form Tr[/(p)] can be com- 
puted through the corresponding function of the eigen- 
values J2i f(^i)- This follows from the fact that the trace 
is independent of the basis in which the density matrix 
is expressed. In the case when the system is entangled, 
knowledge of the spectrum of both the entangled state 
and the reduced subsystems can be used to test for en- 
tanglement through the majorization condition of Nielsen 

lata. 

In all of the above cases the desired measure was ob- 
tained from a single value: the expectation value of the 
swap operator. We now introduce a simple experiment 
to measure (Vat) for a system of N harmonic oscillators 
in an arbitrary state p. 

The experiment, illustrated in Fig. I|T|I. is conducted 
in two stages. First the N oscillators evolve under the 
action of the unitary operator f2. Following that a mea- 
surement is performed in the energy eigenbasis of each 
oscillator and the expectation value of an operator D is 
measured. The specific form of fi and D is such that 

Tr{n P &D}={V N ) (8) 

for all states p. Using the cyclic property of the trace 
this implies 

Vat = Q^DCl. (9) 
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FIG. 1: Apparatus to measure Tr(pVjv). The evolution Q gen- 
erates a discrete Fourier transformation of the -/V oscillators. 
Following that an energy measurement of JV — 1 oscillators is 
performed from which Tr(pVjv) is estimated. 

Since, in general, Vn is not Hermitian, it is not strictly 
possible to associate D in Eqns © and © with an ob- 
servable. Nevertheless, it is possible to express a general 
operator as a weighted sum of POVM elements as was 
done in [l^. From the linearity of the trace, the ex- 
pectation value of the operator corresponds to the same 
weighted sum of measured probabilities. 

From the action of Vjv on the basis state in Eqn. (J5J 
it can be shown that any operator Oj acting solely on 
the j th oscillator transforms as VnO^V^ — > Oj+i (mod 
N). Of particular importance are the creation and anni- 
hilation operators, oj and a,j respectively, which satisfy 

the canonical commutation relation [<Zj,aj] = Sij. Using 
the Kronecker delta function the transformation of the 
creation operators can be written as 

V N a]V^ = J24kj+i- (10) 

i 

Since the set of creation operators {aj} generates an or- 
thonormal basis from the ground state 1 0) <£> 1 0) <£> . . . <B> 1 0) , 
Eqn. I|1U|I . along with the condition 

V N \0) ® |0) <g> . .. ® |0) = |0> ® |0) ® . . . ® |0>, (11) 

is an equivalent definition of Vn- This is a simple but 
important result as it shows that the swap operation gen- 
erates a linear and unitary transformation of the set of 
creation and annihilation operators. We now seek a so- 
lution to Eqn. © where the operators and D are of 
the same kind. Specifically, we require 

QajQ 1 = J2 a i Q » (12) 

i 

Da l Dt = Y,".l } :r (I 3 ) 

i 

where the coefficients fi^ and D^, to be determined 
later, are elements of the unitary matrices and D re- 
spectively. 
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In general, any such unitary operator U which trans- 
forms a set of N creation operators linearly and unitarily 

Ua)tf=Y,4Uii (14) 

i 

is an element of the U(iV) Lie group A property 

of the Lie group is that there corresponds a Hermitian 
operator H of the general form 

ij 

with Ay = A* i5 which generates the group element U — 
exp(-iH). Noting that an arbitrary generator H acting 
on the ground state is identically zero, we can Taylor 
expand U as powers of H and show that Eqn. Ijllfl is 
satisfied for any unitary U of the Lie group, including 
the swap operator. 

The physical importance of Eqn l|12|) is that the evo- 
lution operator Q can be implemented with a linear 
coupling H amiltonian of the form given in Eqn. H15(l . 
iReck et al.l have shown how any such multi-particle cou- 
pling Hamiltonian can be im plem ented using a sequence 
of two particle interactions 20] . In an optics context 
this corresponds to an array of beam-splitters and phase- 
shifters. As for the measurement, our requirement that 
D be an element of U(iV) as- well-as correspond to a mea- 
surement in the energy basis implies that it must be gen- 
erated by an operator of the form Hrj — Y]j ®j a } a j- E x ~ 
plicitly, we require 

D = cxp[— i8ia\ai) ® exp[— «#2<4 a 2] ® ■ ■ ■ 

. . . ® exp[— i0Na/ N a,N]- (16) 

where 0j are free parameters yet to be determined. 

To derive the specific form of the matrix elements rib- 
and the coefficients 0j we substitute Eqn. Q into l|l(J|) 
and derive, with the help of (|12|) and (|13|) . the matrix 
equation 

Vat = f^Dfi (17) 

where [Vjv]ij := 5< j+i. We note that the commutation 
relation between the creation and annihilation operators 
implies 

exp(— itpOjaj) a k exp^aja,-) = a\ exp(— i^Sj^k), (18) 

from which it is straightforward to show from (|16|l and 
(|13|l that the matrix D is diagonal with elements Djj — 
exp(—i8j). The matrix l~i in (|17|) then is such that it 
diagonalises Vjy. Using standard techniques we can solve 
Eqn. (TJJ to give 

n zj = uj^/Vn (19) 

D n = ui~ x (20) 



where to = exp(i2n / N) is the N th root of unity. We see 
that the matrix ft is a discrete Fourier transformation 
and defines the action of the unitary operator Q. The 
unknown phases in (|16[) are found from H20|l to be 

9j = 27r(j - 1)/N, j = l,2,...,N (21) 

and characterize the operator D in (|16fl . To express this 
in terms of POVM elements we rewrite the energy op- 
erator a}-aj in the energy basis of the j th oscillator as 
n j\ n j)( n j\- With this <fR)fl becomes 

D = ^ w i( n i) ■ ■■'W N (nisf)\ni . ..n N ){ni . . .n N \ (22) 
where 

Wj(rij) = exp(—i9jnj) (23) 

is a complex weighting coefficient. In the context 
of a measurement the projector |m . . . njv)(^i • • • nj\r| 
is associated with the joint measurement outcome 
(m, ri2, ■ ■ ■ , tin), which is interpreted as the outcome ri\ 
occurring at the first oscillator, n-i at the second oscil- 
lator and so on. The probability of observing the event 
(ni, ri2, ■ ■ ■ , tin) is given by the overlap of the state with 
the associated projector. For the state VLpStf the joint 
probability is 

Pr(ni, ri2, ■ ■ ■ , tin) — Tr(J7pr2^|ni . . . n.jv)(^i • ■ ■ tin\)- 

(24) 

From lj2"2"fl and lf^I|) the expectation value Tr (Qpfl^D) can 
be expressed as a (complex) weighted sum of measured 
probabilities, 

Tr(n j on + J D) (25) 
= ^2 w 1 (n 1 )w 2 (n 2 ) ■ ■ ■ w N (n N )PT(n 1 ,n 2 , ■ ■ .,n N ) 

which is, by definition, the expectation value (Vat). Inter- 
estingly, 0\ — which means that the weighting function 
Wi(ni) = 1 and is independent of the measurement out- 
come ni of oscillator one. Accordingly, no information 
about the state of oscillator one (after the interaction) 
is used in the estimate of (Vn)- To simplify matters ex- 
perimentally, no measurement need be performed on this 
oscillator. 

To illustrate the practicality of this apparatus we will 
briefly discuss the simplest case which is when N = 2. 
The circuit diagram of this experiment is illustrated in 
Fig. [21 As previously mentioned, for separable states p\ 
and p2 the experiment estimates the overlap Tr(pi/92) 
between two unknown states. By changing the states at 
the input of the apparatus we can choose the quantity 
we wish to measure, ranging from the purity when both 
systems are in identical states pi, through to the fidelity 
(a|p|a) if one of the states is pure \a). On the other 
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p 1 ; 1 50/50 

FIG. 2: Apparatus to measure Tr^pVi). No measurement is 
performed on oscillator one, while detector two measures the 
average value of the parity. The two oscillators interact via a 
50/50 coupling device. 

hand, for entangled states the experiment acts like an 
entanglement witness giving a negative value for some 
entangled states while ensuring a positive value for all 
separable states. 

From Eqns (|12|l and Ijl9(l the required unitary trans- 
formation is of the form 

This specific unitary is generated by an interaction 
Hamiltonian of the form in(a\ci2 — a\a\) applied for a 
time 7r/(4/t), where k is the interaction strength. In op- 
tics, where the two oscillators are realized by different 
spatial modes of the quantized field for example, this 
can be achieved with a simple 50-50 beam-splitter. In 
the case where the oscillators are distinguished by dif- 
ferent polarization and/or temporal modes, then addi- 
tional polarizing beam-splitters and/or time delays would 
be required. Following the interaction a measurement is 
performed in the energy basis of each oscillator and the 
probability Pr(ni,n 2 ) is observed after repeated trials. 
To extract (V2) the distribution is weighted against the 
coefficients w\(ri\) — 1 and u> 2 (n 2 ) = (— l)™ 2 given by 
and lEU- The result is 

(V 2 )=5>irPr(n 2 ) (27) 

where Pr(n 2 ) = J2 ni P r ( n i; n 2), and is independent of 
711, the outcome of the measurement on oscillator one. 
To simplify matters experimentally only the distribution 
Pr + = J2 n2 Pr(n 2 ) and Pr_ = 1— Pr + of the second oscil- 
lator need be measured. The expectation value (V2) can 
then be obtained from the difference Pr + — Pr_, which 
is known as the average parity and corresponds to the 
zero of the Wigner function[2l|. Similar simplifications 
reside in the measurement of (Vn) for higher values of N. 
It is a surprising result that such a range of meaningful 
measures can be obtained from an experiment that re- 
quires no non-linearity, no interferometers, just a linear 
interaction and a single parity measurement. 

We note that higher order moments of the swap oper- 
ator ((V;v) fc ) can also be obtained from the general ap- 
paratus illustrated in Fig^ This can be seen by writing 
(Vat)* 5 = £l Jf D k £l. Physically, this corresponds to the 



same evolution ft of the state p followed by a measure- 
ment of D k which can be written as 

D k = exp[— ik0ia\a\\ <8> exp[— ikQ2a\a,2] ® ■ ■ ■ 

. . . <g> exp[— ik0Na N ciN]- (28) 

Repeating the calculation it is seen that this corre- 
sponds to measuring the same probability distribution 
Pr(ni, n 2 . . . njv), however now the distribution is weight- 
ing by the functions Wj(rij, k) = [wj(rij)] k . 

The measurement procedure introduced here gener- 
alizes to measure the expectation value Tr(pU) of any 
unitary operator U that transforms the creation opera- 
tors of N harmonic oscillators linearly. The key is that 
the associated unitary matrix U can always be diagonal- 
ized as fi^DO where both fl and D are associated with 
the unitary evolution f2 and the measurement D through 
Eqns l|12H and l|13|l respectively. The expectation value 
Tr(pU) is then given by the weighted sum of probabilities 
as in l|25|) . where the specific values of the phases 8j are 
determined from the diagonal matrix D. 

In conclusion, we have introduced a procedure to di- 
rectly measure the quantity Tr(pU) for any unitary op- 
erator generated by a linear coupling Hamiltonian. The 
procedure removes the experimentally challenging step of 
entangling an ancilla to the N systems which was present 
in previous proposals. The result is a more practical pro- 
cedure involving only linear interactions between N oscil- 
lators and local energy measurements. In the case where 
the unitary operator is the JV-particle swap operator, the 
measurement corresponds to a range of useful measures 
depending on the input state, including the purity, the 
fidelity, the overlap, an entanglement witness and gener- 
alized non- linear functionals. 
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